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^ ■ Abstract 

The observed pattern of neutrino mass splittings and mixing angles indicates that their 
family structure is significantly different from that of the charged fermions. We investigate 
the implications of these data for the fermion mass matrices in grand unified theories with 
a type-I seesaw mechanism. We show that, with simple assumptions, naturalness leads to a 
strongly hierarchical Major ana mass matrix for heavy right-handed neutrinos and a partially 
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t^- , cascade form for the Dirac neutrino matrix. We consider various model building scenarios 

OO 

which could alter this conclusion, and discuss their consequences for the construction of a 



natural model. We find that including partially lopsided matrices can aid us in generating 
a satisfying model. 

1 Introduction 

The measurement of neutrino mass splittings and mixing angles [| has provided a new window 
into physics beyond the Standard Model. The fact that the hierarchy between at least one pair 
of the neutrinos is weak and that two leptonic mixing angles are large, in contrast to the strongly 
hierarchical masses of quarks and charged leptons and small CKM mixing, was initially surprising. 
It leads us to surmise that neutrino masses arise through a somewhat different mechanism than 
the quark and charged lepton masses. Thus, the relation between the charged fermion and 
neutrino observables is not necessarily obvious. In fact, we have such a mechanism in the form 
of the type-I seesaw 0, which can naturally yield neutrino masses in the range indicated by 
experiment. Moreover, the physical light neutrino mass matrix is a product of more fundamental 
matrices. This fact can potentially explain the differences between the mixing angles and mass 
hierarchies of the charged fermion and neutrino sectors. 

The seesaw mechanism arises naturally within a grand- unified theory (GUT) such as SO (10) 
H, where each generation of standard model fermions is unified into the 16-dimensional spinor 



1 



representation, together with the right-handed neutrinos. The breaking of B — L (where B and 
L denote baryon and lepton number, respectively), which is a subgroup of SO(10), automatically 
gives rise to Majorana masses for the singlet neutrinos, and thence to the seesaw mechanism. 
Indeed, the neutrino data have encouraged GUT model building || |(| . 

Although GUTs provide a natural framework for massive neutrinos and, combined with family 
symmetries or textures, have allowed for a number of successful models of quark masses and 
mixing, it has proven difficult to incorporate neutrinos in a completely satisfactory manner. In 
this paper, we reconsider neutrino masses and mixings under the guidance of naturalness. That is, 
rather than focusing on a particular theoretical structure and modifying it as necessary to obtain 
the best fit to the data, we will try to minimize the dependence on specific model assumptions 
and work up from the experimental data to see where it naturally leads us. In particular, we will 
show that the construction of a natural, unified picture of all standard model fermion masses 
and mixing angles imposes non-trivial constraints on the structure of both sectors. 

In this framework, we are interested only in the orders of magnitude of various parameters 
and, in pursuing natural solutions, we seek to avoid unnatural cancellations, i.e., that terms 
of a given order must cancel to produce a term of lower order. It may be possible to arrange 
such cancellations in a technically natural way via a judicious choice of symmetries, but this is 
by no means trivial. Furthermore, an exact symmetry is a strong assumption to make, given 
the current uncertainty in the neutrino data. We will instead adopt naturalness as described 
above, seeking to constrain the approximate structure of our theory without ad hoc symmetries. 
Ideally, this structure can serve as a guide for developing well-motivated symmetries upon which 
an ultimately satisfying theory can be built. 

Of course, one must make some assumptions based on previous successes to make progress 
and, in this capacity, we will focus on the SO(10) models with small representations [pi |^]. This 
scenario will serve as a concrete example; however, much of the analysis could be adapted to 
SO(10) models with large representations and/or type-II seesaw mechanisms, as well as to other 
unifying groups. 

This paper is organized as follows: We start by introducing our theoretical framework in 
Section Q and reviewing the experimental data in Section |3|. In Section ^ we derive natural 
constraints on the neutrino mass matrices. Since the fermion mass matrices are related by the 
GUT symmetry, we study the implications of quark mixing in Section [|. In Section |6] we show 
how mass matrices consistent with our constraints can be generated via family symmetries, 
and we investigate how well they can fit the charged fermion masses. In SO(10) models with 
small representations, the neutrino Dirac mass matrix can receive additional contributions via 
couplings to a second up- type Higgs doublet, present in the B — L breaking Higgs field. We 
consider this possibility in Section |7|, supplemented by an Appendix. The remaining sections 
are devoted to two cases which generalize beyond our initial assumptions. These involve models 
wherein otherwise negligible leptonic rotations play an important role in neutrino mixing, either 
due to a lopsided structure in some mass matrices (Section |8|), or to a particular form for the 
effective neutrino matrix (Section [|). We conclude in Section |T(]. 
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2 General Structure of Theory 



The standard model fermions are found in three copies of the spinor representation 16j.[] We will 
make use of the small representations 10//, 45//, 16//, 16' H , 16//, and potentially 16^ to break 
the GUT symmetry and to generate fermion masses. Several authors have used this framework 
to build interesting models J7|, 

The SO (10) symmetry is broken to the Standard Model by GUT scale vacuum expectation 
values (vevs), one in the SU(5) singlet direction of 16// and 16//, denoted v, and (45//) along the 
B — L direction. The electroweak symmetry is broken when weak doublets in 10// acquire vevs. 
It is also possible that the doublets in 16' H and 16'// acquire weak scale vevs, in which case the 
light Higgs doublets are a mixture of weak doublets from the vector and spinor representations 
0. We will assume for now that \Q H does not acquire a weak vev. 

Charged fermion masses are generated via several operators: the renormalizable operator 
16jl6jl0//, which contributes to all Dirac mass matrices for the standard model fermions; the 
higher- dimensional operator 16jl6-,10//45//, which differentiates the quark mass matrices from 
the lepton matrices due to their differing charges under B — L; and 16j 16^16// 16' H , which 
contributes only to down quark and charged lepton mass matrices. The operator 16j 16^10// is 
symmetric in generation space while 16jl6jl0//45// is antisymmetric (16j and 16., are contracted 
as a 120, for (45) oc B — L this is the only contraction that contributes to the mass matrices). 
The operator 16j 16^16// 16' H may be symmetric or asymmetric, depending on how the fields are 
contracted. 

With this set of operators, the Dirac neutrino matrix Md receives contributions from the 
operators 16j 16^10// and 16jl6jl0//45//, and we expect it to be somewhat similar to the up 
quark matrix, i.e., to have a similarly strong hierarchy of mass eigenstates from the first to the 
third generation. For the up quarks this is approximately five orders of magnitude. Although the 
neutrino hierarchy can be somewhat weaker due to factors of 3 coming from the B — L direction 
vev of the 45//, one would still expect roughly a 10~ 4 ratio between the lightest and heaviest 
Dirac matrix eigenvalues. 

We define the orientation of Md as v % M l D ] N\ where N is the Standard Model singlet. Then 
we can parameterize the Dirac matrix as 

M D = L D D D R D . (1) 

Here and throughout the paper the matrices M are dimensionless and the largest eigenvalue 
is normalized to 1. Since we are primarily concerned with interfamily relations this causes no 
problems, but one should bear in mind that there is an overall scale associated with all mass 
matrices. In the above case, the dimensionful Dirac mass operator is u vMpN, where u is the 
mass of the largest eigenvalue. Similarly, throughout the paper L and R will signify unitary 

lr The subscripts i, j will be used to indicate generations while Higgs fields will be denoted with a subscript H . 
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matrices defined by the diagonalization equations 

i:\IM I. = If \ I MR = D 2 = diag {rf , e 2 , l) , (2) 

where rj, e, 1 are the normalized eigenvalues of M. 

In general, and Rn are arbitrary unitary matrices and Dp is a diagonal matrix of the 
eigenvalues of M^; however, we expect the eigenvalues to be strongly hierarchical. This hierarchy 
will be naturally generated if we posit the forms 




D D = diag {v,e,l), 1 p'yfe], R D ~ I pJ% 1 P yfe\. (3) 



We expect ?; < e < 1. Based on the quark hierarchy we may estimate their approximate size as 
rj ~ 10~ 4 and e ~ 1CT 2 , but most of the analysis does not depend on this assumption. 

Ld and Ru are unitary matrices and the parameterizations above should be read as giving 
the orders of magnitude only of the various entries. The parameters p, u, p and their primed 
counterparts are generally expected to be less than or equal to order one. If they were significantly 
larger, various entries would need to cancel to preserve the smaller eigenvalues. Thus p, is, p ~ 1 
is the minimal requirement for naturalness in the absence of an exact symmetry relating the 
Yukawa couplings. This is known as a geometrical hierarchy pattern [ITJ. It corresponds to the 
following form for Mp: 

M D ~ I < e sfe I • (4) 




The central feature of such a matrix is that the off-diagonal entries play a dominant or codominant 
role in determining the two smaller eigenvalues. A geometric hierarchy can be easily obtained 
with a U(l) symmetry via the Froggatt-Nielsen mechanism |TT] . 



On the other hand, p, v, and p may be arbitrarily smaller without endangering the eigenvalue 
hierarchy. In this case the diagonal entries in become dominant and must be correspondingly 
close to the eigenvalues. We will refer to this possibility as a sub-geometric hierarchy. With three 
generations it is, of course, possible to have a mixed case which is partially geometric and partially 
sub-geometric. 

There is one exception to these naturalness considerations, which occurs if is highly 
asymmetric, i.e., if (Mo)ji and (Md)^ are of different orders for some i and j. However, if it 
arises only from 16jl6jl0# and 16il6jl0/f45#, we would not expect this; these operators give 
symmetric and antisymmetric contributions, respectively, which would have to be arranged to 
cancel in a seemingly unnatural way. Thus we generally expect L D and Rd to have similar values 
for their parameters, i.e., p ~ //, v ~ v' and p ~ p' . 

To implement the Type-I seesaw, we need a matrix for the heavy neutrinos: N l M % ^NK Such 
a coupling may arise from — (Mjj)y 16^16^16^16^ when 16// acquires its GUT scale vev v. This 
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non-renormalizable operator is suppressed by some mass m, which is by default the Planck scale 
but which in practice may be somewhat less, depending on the origin of the effective operator. 
The seesaw formula then gives 

M v ~ -M D M~ l (M D ) T . (5) 

As discussed above, M D , M v and M R are dimensionless. The massive parameter which sets the 
scale for the neutrinos is u 2 m/v 2 . For u ~ 100 GeV, v ~ 10 16 GeV, and m ~ mpi ~ 10 18 GeV, 
this comes out to be 0.1 eV, consistent with the range indicated by experiment. 

We stress that the discussion above depends very little on the assumption of small represen- 
tations or the vevs used to do symmetry breaking. One may for example use (45#) proportional 
to the hypercharge generator or use a 54# in place of the 45# to accomplish the breaking from 



SU(5) to the standard model |jl2fl . Alternatively, we could have used the large representation 
approach with 10#, 120//, and 126#, which many authors have used for model building [ |T3|j . 
In any case, we still expect a hierarchy in the quark and charged lepton mass matrices. Due to 
SO (10) relations, this hierarchy should manifest itself in the Dirac neutrino matrix as well and 
the same naturalness considerations apply. 



3 Experimental Constraints 

The detection of neutrino oscillation is successfully explained by massive neutrinos with non- 
trivial mixing. We know two mass squared splittings among the neutrinos and two mixing angles 
of the leptonic mixing matrix, with a limit on the third for the physical light neutrinos [[J, 

tan 2 12 = 0.45 ± 0.05 ; Am s 2 ol = (8.0 ± 0.3) x 10~ 5 eV 2 ; 

sin 2 20 2 3 - 1.02 ± 0.04 ; A™at m = ( 2 -5 ± 0.2) x 10~ 3 eV 2 ; 

sin 2 26 13 = ± 0.05 . (6) 



Additionally, cosmological considerations place a limit on the total mass of the neutrinos 
along with limits from tritium beta decay and neutrinoless double beta decay on the electron 
neutrino |l|, |2], 15]. These experimental results constrain the total mass of the light neutrinos to 



be less than or of the order of 1 eV. Our discussion does not depend on the exact number since 
the masses are degenerate in this limit. The bound will only become important to our analysis 
if it approaches the atmospheric mass splitting. 

The mixing is characterized by the PMNS matrix, a unitary matrix parameterized by three 
angles and three phases, 

Vpmns = L\L U (7) 

C12C13 S12C13 s 13 e~ lSs 

-S12C23 - c 12 s 23 s 13 e iS C12C23 - s 12 s 23 s l3 e iS s 23 c 13 ] x diag (e tai/2 , e ia2/2 , l) . 

S12S23 — Ci2C23Si3e l<5 — C12S23 — Si2C23Si3e J<5 C23C13 
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For concreteness, we will assume the tribimaximal solution which sets the mixing angles 
0i 2 = arcsin (l/Vs) ~ 35°, 13 = 0°, 23 = 45° 0, 



PMNS 



V 



\ 



V 



(8) 



neglecting phases. This is in some sense an extreme solution consistent with the data. Given 
the several seemingly disparate factors which influence the angles, it seems highly unlikely that 
any model will predict exactly zero for 6 13 , or exactly maximal atmospheric mixing, unless 
carefully designed to do so ||. Therefore it may well be that experiments eventually favor a less 
striking set of angles. Furthermore, in a detailed model one would also need to carefully consider 
renormalization, which can have a significant effect on the mixing angles and mass splittings [ F7 
We do not address these effects in further detail in this paper because they make little difference 
in our analysis. We are only looking at relative orders of magnitude of masses and mixing angles. 
Due to its simple structure, we will use the tribimaximal solution as an experimental input. The 
critical facts we need are the existence of two large neutrino mixing angles and a relatively weak 
neutrino mass hierarchy, both of which will remain true despite renormalization effects. 

We will assume for now that the tribimaximal structure is generated essentially in the neutrino 
sector; given the charged lepton hierarchy, we usually expect relatively small rotations in L e 
compared to the large PMNS entries. Since we are only concerned with orders of magnitude, we 
will (for now) neglect the charged lepton component. As with the geometric hierarchy discussed 
in Section [| there is one exception to this rule associated with a highly asymmetric structure, 
this time in the charged lepton matrix. Such a lopsided matrix can introduce large rotations, as 
shown in the Albright-Barr model ]/]]. This case will be discussed further in Section 

The neutrino mass matrix will be diagonalized by the tribimaximal rotations if it has the 
form 



VpMNS^fVpMNS OC 



{nil + \ m -i) 
-| (mi - m 2 ) 
~ (mi - m 2 ) 



- 5 (mi - m 2 ) 
~ (\mi + m 2 + |m 3 ) 
-§ (~mi + m 2 - |m 3 ) 



\ (mi - m 2 ) 
-| (|mi + m 2 - |m 3 ) 
\ (|mi + m 2 + |m 3 ) 



(9) 



i.e., L v = R v = Vpmns- The m's are the physical neutrino masses with an arbitrary phase for 
mi and m 2 . Since we know the two mass squared differences, we may rewrite these in terms of 
a single mass, 



m x 



e^ 1 |mi| 



m 2 



e^Vlmil 2 



A 2 

^sol ) 



m 3 



mi| 2 + A s 2 ol ±A 2 tm , 



(10) 



2 For example, a bimaximal mixing scenario (6*12, 6*23 = 45°, 6*13 = 0) at the GUT scale can produce weak scale 
mixing angles consistent with the data quoted above |l7| . 
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where we have introduced the notation A = V Am 2 . The ± in the definition of 777.3 represents the 
choice of normal (+) or inverted (— ) hierarchy. We take the phase factors e 1 ^ 1 ' 2 to be ±1 so that 
there are just a few choices of relative positive or negative to make. Since we are only concerned 
with orders of magnitude and this will give the extrema, this should not limit the analysis. Then 
it is simple to scan through the allowed range of mi. By doing this, one can observe the patterns 
of relative order in the neutrino entries which are consistent with experiment. The potentially 
interesting possibilities are 

1. M v ~ ||j 1 Y corresponding to mi <C m 2 ~ A so i, normal hierarchy. 

2. M v ~ fSjjj, corresponding to 2mi ~ m 2 — ^A so i, 2 — 4>x = t, normal hierarchy. 

3. M u r>j (011), corresponding to A sol (A atm ) < m 1 ~ m 2 < A atm ( v / 2A atm ), <p 2 - ^ = 0, 
normal (inverted) hierarchy. 

4. M v ~ ?)' corres P on ^ n S ^° degenerate masses, 2 = 0, 0i = 0. 

5. My ~ f 1 Y corresponding to degenerate masses, 2 = 7r, (j>i — n. 

6. My ~ (|]])' corresponding to degenerate masses, </> 2 — <pi = ir. 

Here A = ~ 0.2 and should be read as at least a few orders of magnitude smaller than 1. 
Any other possibilities should be roughly an interpolation between those listed and we do not 
expect them to lead to significant deviations from the results following. 

The cases with non-degenerate masses, namely the first through third above, violate the 
geometrical hierarchy naturalness limit discussed in Section |2|. In each case the democratic 2-3 
block generically leads to two large eigenvalues of order 1 and one large mixing angle. Then the 
couplings of the first generation give a naive estimate for the third eigenvalue of A, A 2 , and 1 for 
the first, second, and third cases, respectively. This is not compatible with the eigenvalue ranges 
listed above, so some unexpected cancellations would have to take place. Moreover, these cases 
are more compatible with a small 812 due to the smallness of all off-diagonal first generation 
entries. The fourth and fifth cases naturally lead to degenerate eigenvalues as listed but imply 
unnatural precision to account for the large mixing angles. 

In short, hierarchical neutrino masses are unexpected in conjunction with large mixing angles, 
and large mixing angles naturally proceed from large off-diagonal entries in the effective mass 
matrix. Thus, case 6 above is the most natural simple assumption to account for the experimental 
data; it is known democratic mass matrix JT8| . 



We can also consider evidence from neutrinoless double beta decay experiments. A positive 
signal would confirm the Majorana nature of neutrinos and lend credence to seesaw models. The 
experimental status is controversial: After the Heidelberg-Moscow collaboration set the limit 
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\m ee \ = \(M u ) n \ < 0.35 /i eV, where h denotes the uncertainty of the nuclear matrix element 
P, [151, a subset of the collaboration claimed evidence for a signal [|19[]. Depending on the value of 



h, this signal points at quasi-degenerate neutrino masses in the range 0.1 — 0.9 eV [[|. This result 
clearly requires confirmation from current and future experiments. If confirmed, the hierarchical 
scenarios would be ruled out, consistent with our conclusions from naturalness. However, since 
this claim is still controversial [E0| , we will not rule out the hierarchical scenarios in our analysis. 



We note that for an inverted hierarchy with ni2 — ^=A atm , we could have 

/I 1 IN /ll IX 
M v ~ 1 11, 110, (11) 





depending on the phases 0i 5 2- These should be thought of as special subcases of case 6. As will 
be shown in the next section, these possibilities will only add additional modeling constraints 
compared to case 6 without additional explanatory power, so they are not particularly interesting 
in this context. Bearing these caveats in mind we shall, however, consider some cases besides 6 
because they may relax other naturalness constraints. 



4 Modeling 

Now we will do a little rearranging of the seesaw formula in terms of the eigenvalues and unitary 
matrix decomposition of Mr,: 



D-Jl} D M v U D D- D \ 



(12) 



Applying this to the sixth and henceforth canonical case above, we get 



/_L L 

„9 



D 



i/c 
1 



1 



(13) 
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where we have kept only the leading terms. The salient point is that, with the assumption 
/A v ' i p' < 1, the Ld rotations (and similarly the charged lepton rotations) cannot change the 
orders of the entries. From this we see the apparent double hierarchy for Mr. its eigenvalues 
naturally scale as rf, e 2 , 1 compared to r], e, 1 for Mo- 

Most of the other cases are similar and retain at least a -4 ratio between the first and third 
eigenvalues. For the cases where M v has entries less than order one, the unitary rotations can 
contribute significantly, in particular they can "fill in" the zero entries, but they cannot make 
any entries larger than order unity in 
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There are two cases which may differ importantly from the others. Case 1 in Section |3| is 
interesting since it yields 



R ] D M- l W D ~ 



Similarly, for the second case we get 



V 



i(A + A *'fl) i(A + ^)\ 
i i 

i 1 / 



(14) 



(15) 



In these cases we see that we have mitigated the largest ratio of entries from 4 to a smaller 
value, although said ratio remains significantly larger than 

Let us now consider the effects of the matrix Rp on the canonical case. We will show that, 
under the current assumptions, one can put additional constraints on p, v and p. To begin, we 
parameterize the inverse heavy neutrino matrix 




(16) 



and evaluate both Eq. ( |P2| ) and 



Mn 1 



R D Do X l} D M v V D Dn X R T D , 



(17) 



which is just another rearrangement of the seesaw formula. Keeping only potentially leading 
terms, we find 



.4 
B 
C 
D 
E 
F 



1 



3/2-rl/2 



+ 



V 



rje 



3/2 



+ 



rje 



1/2 



n 



Tje 7^1/2g3/2 



\1V 1 



c3/2 



V 



+ 



2/i 



2pu 1 



,1/2 



\tv v + pp 



pise 



v 2 

— + 

rj rj 



1/2.1/2 



2pz/ 

1/2.1/2 



+ 



IT/2 



e 



X / 2 p 1 

+ — h - 

1/2 ^ e 3/2 ^ £ ' 



rl/2 



+ 1 . 



;is) 
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Now, with a little consideration, one can see that each entry should only be as big as the rightmost 
term. This is because Eq. (0) must still be satisfied looking only at the order of the terms. For 
example, we can look at the equation for the (12) entry of Eq. (O) in terms of A through F and 
fi, u, and p via Eqs. (0) and flTSI). This comes out to be 



A + B + C (^p^fe + f*v^=) + D + E Op + v) Vv + F P"VV£ ~ ^- ( 19 ) 

B appears in this equation with a coefficient of order 1, thus any solution to the set of conditions 
in Eqs. ( |18D with B > — will apparently not satisfy Eq. fll9l).Pl This is a naturalness condition. 
One can, of course, numerically satisfy both equations but it requires a cancellation between two 
terms to at least an order of magnitude. If we want to avoid the need for a symmetry precisely 
relating various parameters, the only natural solution is to set B ~ — .[] 

Applying the same analysis to the rest of Eqs. (|TB1), we come to the conclusion that 



M^ 1 ~ D D X L D M V L D D D \ (20) 

or that the hierarchy of M^ 1 could be even stronger, regardless of Rd- Then we must impose 
constraints on the mixing parameters in Eqs. (|1J|) so that the parameters B — F do not become 
too large: 

P<yf, "<Vv, P<v^- (21) 

For r) ~ 10~ 4 and e ~ 10~ 2 , this corresponds to fi, p < lO^ 1 and v < 10~ 2 . 

If we take the minimum required suppression and apply it to //, v\ and p' as well, we get the 
cascade hierarchy pattern [TO, |21j for the Dirac matrix, 




(22) 

\T] e 1/ 

For any hierarchical texture of R^M^ 1 R* D we will find that M^ 1 generally retains the same 
hierarchy. Intuitively, this is because Rd will tend to smear out any hierarchy in M^ 1 ; the 
larger entries will be rotated into the smaller. The hierarchy would only be sharpened if there 
were a very precise relation between Rd and M^ 1 , which we have no reason to expect. So in 
general, if RpM^ 1 R* D has a hierarchy of entries, M^ 1 should have at least as strong a hierarchy. 
Conversely, to maintain a strong hierarchy in R D M^ l R* D ^ the unitary rotations cannot be too 
far from diagonal, a fact reflected in the constraints on /i, v and p. 



3 Here the important number is actually the ratio B/F ~ 1/(775). Using the conventions above we find F ~ 1, 
but there is an overall numerical factor which we omit because it can be absorbed into the dimensionful vevs. 

4 Technically, it could be smaller since Eq. ( |l^ ) depends on experimental numbers. Thus in Eq. (|l^), it 
may cancel the theoretical parameter term ]i,j(jf'l' 1 yfe) without fine tuning as long as it is consistent with the 
experimentally allowed range. At any rate, it would only make the hierarchy stronger since A ~ -K regardless. 
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For the other possible textures of M v with one or two suppressed entries, we mostly find equal 
or stronger constraints on the p, v, and p. For example, if the (23) and (32) entries of M u are 
small so that the corresponding entries in l) D M v L* D D~^ are much less than -, then we also 
require E <C -. This in turn imposes stronger constraints on the mixing parameters. This is the 
situation for cases 3-5 as well as the special sub-cases of 6 mentioned in Section [3|. 

It is interesting that the constraints on p, u, and p remain valid even if we take the first case 
of the list, 




/4 

'n 

A_ 

■qe 



■i/c 
1 
*2 



n 1 

i 



(23) 
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This is because we retain the strong hierarchy along the first column and row, as well as in the 
(23)-block, whose entries remain less than or equal in order to the first generation entries. 

The one exceptional case is the other form noted before, case 2. This leads one to the 
conclusion 



Ml 



V J(a+mV?) 



i 

e 2 
1 
c 



(24) 



/ 



and the naturalness conditions 



< 



< 



A 



0.1, 



p < \Te ~ 0.1 



(25) 



So in this case we are not as constrained as the cascade pattern but still more constrained than 
the geometric pattern; only the constraint on p remains the same. This makes sense since, in this 
case, we have a relatively weak hierarchy in the first row and column compared to the canonical 
case. Therefore, we find weaker constraints on the rotation parameters for the first generation. 

In general then, we are led to both a double (or at least enhanced) hierarchy for Mr and a 
cascade (or sub-geometrical) pattern for Mjy in a simple type-I scenario. Other authors have come 
to similar conclusions following from the assumption of hierarchical Yukawa matrices [p2, ^3 



5 CKM Constraints 

The Dirac mass matrices of quarks and leptons are related by SO(10) and possibly family sym- 
metries. Thus, we should also consider the size of the unitary rotations in the up and down quark 
mass matrices, which are measurable through the CKM matrix, Vckm = L^L^. The experimental 
CKM values are Q 

/ 1 0.226 ±0.002 [4.3 ± 0.3] x 10~ 3 \ 

Vckm = 0.23 ± 0.01 1 [4.2 ± 0.06] x 10~ 2 . (26) 

\[7.4±0.8] x 10~ 3 3.5 x 10~ 2 1 / 
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If we suppose for the moment a geometric pattern for both the up and down quark matrices, 
then the predicted CKM matrix is 



V< 



CKM 



V 



+ 



ll'd 

m s 

I m d I / " L u 

m b \l mt 



m c m d 
m t m b 

m c m d 
mtm s 



i rn d I rriu_ _|_ 



' EM _|_ / lii£ 
m b \/ mt 



m u m s 
m t m b 



m u m d 
mtm a 



( 1 

-0.23 + 0.06 -0.001 
V-0.03 + 0.003 - 0.01 



0.23 - 0.06 + 4 x 10~ 4 
1 

-0.14 + 0.04-7 x 10~ 4 



0.03- 
0.14 



m u m s 
in, \l m c m b 

m u m d 
m c m b 



+ 



0.003 - 0.008 N 
- 0.04 + 0.002 



(27) 



nicely reproduces the 



A few features are striking. One is that the geometric ratio 

experimental value for the first-second generation mixing The dominant term comes from 
the down quark mixing, while the contribution from the up quark mixing is significantly too 
small to account for the mixing by itself. Secondly, the down quark contribution to the first- 
third mixing is too large by roughly an order of magnitude. Lastly, the down quark contribution 
to the second-third generation mixing is also too large by roughly a factor of three. So the 
geometric hierarchy does a good job for the Cabibbo angle but gives too much mixing with the 
third generation. 

This result is consistent with a partially cascade structure in Md and M U .Q The relatively 
large Cabibbo angle indicates that the down quark matrix should be close to geometrical in the 
1-2 block. However, it will fit the data better if it is cascade-like in the third generation. If the 
same were true of Md, we would be consistent with the second case from Section [3| On the other 
hand, since 16jl6jl6#16# only contributes to the down quark and charged lepton matrices, the 
neutrino matrix could remain completely cascade-like without conflict. 



6 Implementing the Cascade Hierarchy 

Since we argue that a cascade texture is theoretically desirable, we will investigate how it can 



be generated. We will make use of the Froggatt-Nielsen mechanism ]TTJ and consider a global 
U(l) x7j 2 x Zj' 2 symmetry. We introduce three SO(10) singlets 0j. The flavor symmetry is broken 
spontaneously at a high scale m by vevs of the singlet fields, which we expect to be all of the 
same order, (</>). The symmetry breaking is assumed to be transmitted to quarks and leptons 
through interactions with heavy particles so that the Yukawa couplings are constructed out of 
powers of C = (<p) /m with a texture dictated by the family symmetry. 



5 Since the largest terms come from the down quark sector, the CKM values are also consistent with a geometric 
hierarchy in M u . Given SO(10) relations and possible family symmetries, the simplest assumption is that M u has 
a similar hierarchy structure to Md- 
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We assign the following charges: 



Field 
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16 2 
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z 2 
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Then the operator M[q 16^16 jlOn originates from ^16il6jlOH, where $ represents the higher- 
dimensional couplings, 



in 



'2 



■ (0l) 3 03 
h (0l) 2 
X 201 02 



>2 <?3 



■m- 



^2 0102 
1 



(2f 



so that 



M- 



10 




(29) 



This is the cascade form of Eq. ( p2|) with n = ( 4 and e = ( 2 . The same pattern can easily be 
reproduced in the other operators which contribute to fermion masses. Note that in the absence 
of the 7j2 symmetries we would have generated a geometric hierarchy. 

We must also consider whether a cascade hierarchy can naturally accommodate the fermion 
masses in a unified theory. Restricting ourselves to two generations, the operators discussed in 
Section |2] contribute to the (normalized) mass matrices as follows: 



a' a + (3 
a- (3 1 

a' + i a + + 7 
ct-/3 + 7 1 



M D 

M P 



a 1 a — 3(3 
a + 3(3 1 

a' + i a - 3(3 + 7 
a + 3(3 + 7 1 



(30) 



Here, the terms a and a 1 parameterize the operator 16jl6jl0#. The parameter (3 derives from 
16jl6j 10#45#, while 7 and 7' characterize 16il6jl6 hIG'h- Looking at the determinants, we 
calculate the mass ratios: 



m 



c I / 

~ a 



7Bt 



a 2 + /3 2 



a: 



a 2 + 9(3 2 



— ~ la' + y + Z? 2 - (a + 7)' 



n? 



— t ~ U' + 7 ' + 9^ 2 _ (a + 7 )1 . 



(31) 



As expected, /5 accounts for the difference of down quark and charged fermion masses, 



8P 2 



m,. 



4 x 10" 2 - 



x 10 



-2 



+ 



(32) 
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where we used (m /i /m r ) GUT ~ 0.06 and (m s /mf,) GUT ~ 0.02. Since we wish to minimize off- 
diagonal terms in a cascade-like matrix, we will use the smaller value for /?,[] 

P ~ 7 x 10" 2 . (33) 

Then we obtain 

e = ^ + 8/ 3 2 ~7x 10- 2 , (34) 

with (m c /m() GUT ~ 0.03. 

In order to have a cascade form for Mc, we require a' ~ a ± 3/? ~ e. Since 3/? ~ 0.2, this 
implies a' ~ 0.1, independent of a. This value of a' can be consistent with the value of e in 
Eq. (|34"D, but it needs to cancel significantly with a 2 to ensure a suitably small value for m c /m t . 
Conversely, m c /mt implies a' < 10 -2 , which leads to a geometric hierarchy in Mp. Since we 
have been trying to avoid requiring the cancellation of theoretical parameters, this simple cascade 
ansatz is problematic. 

One particularly attractive way out of this dilemma is to consider the possibility that Mp, 
but not M u , receives additional contributions, e.g., via particular higher- dimensional operators. 
If such an operator gave a contribution to the (22)-element of of order e ~ 0.1, a' could be 
made sufficiently small. We consider such a scenario in the following section. 

7 New Contributions to Mp 

In Section |] we saw that the observed pattern of neutrino masses and mixings leads us to an 
enhanced hierarchy for Mr, compared to Mp. One should note, however, that while Mp is related 
to the observed quark and charged lepton hierarchies by SO(10) and any family symmetries, it 
is not directly observed. In particular, one may include another operator, 16il6jl6Hlf>H- As 
noted above, the weak doublet in 16 H can acquire a weak scale vev v! such that this operator 
potentially contributes to the up quark and neutrino masses. However, it can be constructed to 
contribute only to the Dirac neutrino matrix. In this case we expect u' < u, since u is required 
to generate a large top quark mass and the sum of the squares of weak scale vevs must equal 
(246 GeV) 2 . 

A simple possibility for generating this operator is to integrate out SO(10) singlets, S, at 
some scale above the relevant GUT scale vevs. For this purpose we can propose the operators 

MtfieiieffSj + M^-ieae^ + m s (Ms)^ sSj ■ (35) 

We assume at least three singlets to guarantee that all three righthanded neutrinos become heavy. 
As usual, we define m s to have units of mass so that Ms is dimensionless with entries of order 
1 or smaller, and similarly we normalize M and M in Eq. (^). In the following analysis we 

6 The larger value, (3 ~ 0.1, leads to e ~ 0.1. 
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assume that all the S singlets are integrated out to generate an effective Majorana mass for the 
iV's. To compute this via a straightforward seesaw mechanism, we will work in the basis where 
Ms is diagonal and impose the conditions 



(Ms)u > — 
m s 

for all i. 

The mass matrix for the electrically neutral particles readsQ 



(36) 



{y N S) 
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\uM T D 
\\u'M" T 



{uM D 


\vM T 



M'\ 



\u'_ 

\vM 



rn 



Ms J 




(37) 



As derived in Appendix the light neutrino mass matrix is then given by 



My ~ 



M D [M 



M S M l Ml - | 



MM ^Ml + M D [MM 



T 



u'v 



X 



u m s 



< 1 . (3* 



The mass of the heaviest neutrino is of order u 2 m s /v 2 . It is crucial that, in the final formula, 
M D appears in all terms, i.e., terms quadratic in M M have not appeared. 

Let us study the effect of the new contributions. We parameterize the various matrices as 
follows: 

'A B C> 
B D E 
,C E F i 



M 



M~ 



M S M~ 



m'm 1 
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b 




1/ 
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el 


I 



(39) 



(note that the matrix MM is generally not symmetric), and 



M,/M M 



+ 



m'm 1 



A' 


B' 


a 


B' 


D' 


E' 


a 


E' 


F' 



(40) 



The last matrix, with primed capital letters, is the total effective matrix which takes the place of 
Mft 1 in Section |j. The unprimed capital letters parameterize the familiar heavy neutrino matrix 
and the lower case letters parameterize the new terms. Before proceeding to consider the effects 
of these new terms, we note that Mr can easily acquire a double hierarchy if it is generated by 
integrating out heavy singlets, as described above. If M has a hierarchy comparable to Md and 
Ms is roughly democratic, a double hierarchy occurs naturally. 

We can write the total effective parameters in terms of these old and new components and 
perform the same analysis on the total effective matrix [A 1 — F') as we did on the simple type-I 
parameters (A — F) in Section ^. Then we obtain the following set of equations: 

v' 

_7] 



A' ~ A + 



+ c r 
rje 



x 



7 Barr calls this scenario a type-Ill seesaw mechanism |25); however, it can also be understood as a product of 
two type-I mechanisms. 
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. x/V Ve 
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x 



x 



(41) 



In these equations we have kept only the leading terms. In doing so, we make use of the 
important fact that the constraints on /i, v, and p still apply. They follow from consideration of 
the experimental data and the geometric constraints on Mo only.0 

Although these equations still appear somewhat complicated, the requirement that we fit the 
same hierarchy of orders as imposed in Eqs. (|18j) can only be satisfied in a few ways. In general, 
the new terms give us new parameters which could play a role in a precision fit to the data, but 
they will not affect the conclusions of this paper unless they dominate over the old terms. Let 
us consider the canonical case, which implies 



'A! B' C'^ 
B' D' E' 
,C E' F'j 



oc 



(\ 

V 

J_ 



rje n 

1 1 

? 7 
l 



(42) 
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Examining Eq. (P4|), this puts some initial constraints on our new parameters. For example, 

bx , F' 

— < B' . 43 

i] rje 



These constraints may be summarized in matrix form: 



/ a 


b 






A 

'/ 
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b' 
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e 
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>i 
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ef 


si 






1 
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V 



F' 

x 



(44) 



Taking these restrictions into account, we conclude that to satisfy ~ \ we must have 



.4 



r/ 2 



or a ~ F 



rjx 



(45) 



The latter case is initially appealing because one can apparently trade the strong double hierarchy 
constraint on Mr for a weaker standard hierarchy in M M if the term involving a dominates. 



8 This would not be the case if there were new terms in the effective total matrix which did not involve Md- 
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This turns out not to be feasible. Recall that M and M have all entries of order 1 or less. 



Thus, if a = F— , there exists some % and some n > 1 for which 

, l i n F' 

M u = - and M a = 

n 7]x 



(46) 



(cf. Eq. (^)). Then the assumption that the a term dominates over A gives us the inequality 



1 ax A 

rv, > 



7/- 



T]F> 
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~F' 
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TJX 



{M S ) U F' 



(47) 



2 

from which we obtain (Ms)uF' < Applying the seesaw constraint on Ms and inserting the 
definition of x gives us 



m.. 



-F' < (Ms)uF' < 



u' 2 v 2 
n 2 u 2 m 2 s 



This requires F' to be too small, that is, 



1 < 



M M S M~ 



F < F' < 



J 33 



11 
U V 

n 2 u 2 me 



(Ah 



(49) 



Since v -C m s and v! < u, this condition cannot be satisfied. Thus the additional contributions 
cannot dominate over the type-I contributions or change the need for a double hierarchy. 

One can instead look at case 1 from Section [| If the new terms dominate in the largest ratio, 
which is still A'/F', this implies M a = A 1 ^-. Proceeding as in the canonical case above, one 
finds X(Ms)uF' < ^4. Since we require F' > 1 this is only possible if 



x 2 > (M s )ii\n 2 > — \n 2 



or equivalently, 



A < 



u' 2 v 
n 2 u 2 m. 



< 1. 



(50) 



(51) 



This is a very marginal case since we are relying onu < m s to use the seesaw formula as a valid 
approximation and A ~ 0.2. 

If we proceed nonetheless, then we impose the conditions on B'\ 

X 

rje 



B B 
— > — ~ 

F> ~ F> 



(52) 



which implies the constraint M i2 < xj [ne(Ms)u]- Now we turn to D' ~ -j. By similar reasoning 
as in the canonical case it can be shown that D must dominate to satisfy D' of the appropriate 
magnitude, due to the suppression of the new terms by x. Then 



> 

F' ~ F 



P- 



(53) 
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for some j ^ i, which gives us the condition 



This in turn implies 




(55) 



so as not to violate the bound on b. We find then that the new contributions can technically 
dominate in the (11) entry but the type-I terms remain comparable and dominate in other entries, 
still exhibiting a strong hierarchy compared to Mjj. 

The related case 2, with (M u )n <C A, is, not surprisingly, similar. One finds that the new 
term a can dominate if A A /^ < jO^ , which provides somewhat more room for consistency with 
the seesaw approximation. The constraints on the matrices are 

In both cases the new terms can dominate in some entries, but the type-I terms remain 
important and retain a strong, albeit not quite double, hierarchy. We note that this is due largely 
to the structure of the theory: if Mr is a dimension-five operator generated by integrating out 
singlets, then a hierarchy in M similar to that in Mo naturally leads to a doubled hierarchy in 
Mr. Due to the suppression of the new terms by v/m s , Mr will always play an important role. 

It is interesting that even with the 16j 16^16// 16 H operator only contributing to the neutrino 
sector, we still derive the cascade constraints. Although this operator only contributes to the 
Dirac neutrino matrix, the constraints apply to the operators which generate the up quark matrix. 
This follows from the precise relations between the higher dimensional operators induced by their 
common origin. These relations result in Mo appearing in all terms of the formula for M v . As 
a consequence of the persistent cascade constraints, we cannot use the new terms to solve the 
mass splitting problems discussed in Section |[ 

If one treats 16jl6jl6//16 H and 16,; 16j 16^16// as independent, it is possible to relax said 
constraints. That is, in the discussion above both operators depend on the coupling Mij 16il6#S./ 
and are therefore related. If we allow them to vary arbitrarily, then the modified seesaw formula 
in Eq. ( |38| ) would have additional terms which did not involve Md- In effect, we would be adding 
new terms to the Dirac neutrino matrix which could strongly alter its hierarchy compared to the 
quarks and charged leptons. If this resulted in a relatively weak Dirac neutrino hierarchy, Mr 
would have a correspondingly weakened hierarchy and the mixing parameter constraints would 
also weaken. However, as shown above, this is not necessarily the case when one begins with a 
more complete theory. 



18 



In general, if one can weaken the Dirac neutrino hierarchy without upsetting the charged 
fermion hierarchies, the requirement of a double hierarchy in Mr and a cascade hierarchy in 
Mo becomes less restrictive, since they are specified relative to the eigenvalue hierarchy of Mp. 
One possibility for doing so may be to introduce a vector-like fourth generation of down quarks 
and leptons at the GUT scale. This can relate Mp to the down quark hierarchy such that the 
hierarchy of Mr is similar to that of the up quarks | |26| . 

8 Lopsided Models 

Thus far, we have not allowed for any cancellations between terms in our equations, in keeping 
with our aim to eliminate unnatural models. There are, however, two scenarios where one must 
be more careful. These are cases where unitary rotations play a very significant role either due 
to large rotations or small entries in the neutrino matrix. 

In this section we will consider the first type of these cases, lopsided models, wherein the 
operator 16jl6j 16^16'^ is constructed so as to contribute in a highly asymmetrical way to the 
down quark and charged lepton mass matrices [|7| .[] These lopsided matrices can yield a natural 
hierarchy while violating the geometric pattern limit discussed above. Lopsidedness results in 
large off-diagonal terms in the unitary rotations on one side of the matrix but not both. 

To illustrate these features we will restrict ourselves to two generations first. The follow- 
ing table summarizes the three natural cases we have discussed for a generic matrix M with 
eigenvalues e and 1, which is diagonalized by the unitary rotation matrices L and R. 



Hierarchy 
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(el) 
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(el) 


Lopsided 


(l l) 


(ID 


(ID 



For both the geometric and cascade cases L and R are similar to each other. As expected, the 
off-diagonal entries of L and R for the cascade case are smaller than in the geometric case. 

The lopsided case, being highly asymmetric, leads to very different rotation matrices on the 
left and right. We see that to generate large mixing on one side, i.e., R with all entries of the 
same order, we are led to L being closer to diagonal than in the geometric case. Rather, it is 
similar to the cascade rotation matrices. So in this simple case, to preserve naturalness, there is 
a tradeoff between the left and right sides. If one side's unitary rotation violates the geometric 
naturalness bound, the other's is concomitantly constrained to be closer to unity. 

To take potentially large mixing in the charged lepton sector into account, we have to reeval- 
uate our seesaw formula. In Eq. (|12j), we neglected the rotations from the charged lepton sector, 

9 We will not discuss the origin of these lopsided matrices, which, e.g., can be due to family symmetries @. 
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parameterized by L e (cf. Eq. ([?])). To include them we rewrite the formula as 

R ] D M~ l W D = D^VoM'XDd, 

where 

V = T] 



M' 



PMNS-D^^PMNS 



L D L e 



(57) 



(5f 



M' v is the light neutrino mass matrix in the basis where the charged leptons are diagonal. It can 
have the same forms as discussed in Section || for M v . With the substitutions M v — > M' v and 
L D — > V , the equations used above are unaltered. 

The crucial difference is that the assumed form of in Eq. (|3|) does not necessarily apply 
to Vo in the lopsided case. Since Vo contains off-diagonal entries of order one, we may arrange 
for terms of equal order to cancel each other in VqM' u Vq '. This is not fine tuning because 
we are, in effect, canceling an experimental term with a theoretical one, rather than canceling 
two theoretical parameters against each other. To put it another way, we are simply using a 
theoretical term to generate an experimental parameter of the same order. The result is that 
we may be able to have a form for R^M^ 1 R* D which does not have such a strong hierarchy, and 
which in turn may not imply the restrictive cascade form for Mjj. In such a scenario, some or 
all of the large mixing in Vpmns comes from charged lepton unitary rotations. 

To examine the lopsided case further we must see what can be said about the matrix Vq. 
Again, we can look at the CKM matrix for possible constraints. It can tell us about the potential 
lopsidedness in the down quark and charged lepton mass matrices. The operator we are using to 
generate lopsidedness contributes to M e as the transpose of its contribution to M d , as is familiar 
from SU(5) models.^ Hence, large rotations in L e would coincide with large rotations in R d and 
vice versa. We now see that the experimental values are consistent with either a cascade structure 
or a lopsided structure for the down quark mass matrix in the third generation couplings. 

One might hope that the relatively large 1-2 mixing, which is consistent with a geometric 
hierarchy in the down quark matrix (cf. Section |5|), would constrain the 1-2 mixing in R d . This, 
however, turns out not to be the case. We can construct a matrix with all the desired features 
and generically large righthanded mixing, e.g., 
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(59) 



Thus, although the CKM matrix is highly suggestive of either a partially cascade or lopsided form 
for the down quark mass matrix, it is difficult to constrain the form of R d and its counterpart 
L P in the latter case. 



10 This is simply due to the fact that W H breaks SO(10) to SU(5), so 16^6^16^16'^ is basically an SU(5) 
Yukawa operator for down quarks and charged fermions, suppressed by v/M. 
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On the other hand, we note that since Vpmns has a small value for the (13) entry, naturalness 
requires that at least one entry in the column L* 1 be correspondingly small. This suggests that 
we can rule out the extreme lopsided case shown above. 

Lopsidedness also modifies the eigenvalue fitting we did in Section |6|. Let us consider the case 
where 16jl6jl6#16' H is lopsided and assume that it contributes to only one off-diagonal entry 
in Mrf and M e in a significant way. Then the mass matrices in Eq. are modified to 



M r , 



a' + i a + (3 
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with the corresponding eigenvalues 
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This yields 
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(60) 
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(62) 



I + 7 2 

which is only a slight improvement over the symmetric case, a cascade structure in Mr, is still 
inconsistent with the charged fermion hierarchies. Thus, in the absence of additional contribu- 
tions to the mass matrices, it seems we must rely on large, lopsided mixing between the second 
and third generations to alleviate the need for a cascade structure in the 2-3 block of Mp. 



9 Small Entries and Mixing 



Aside from lopsided matrices, there is another scenario in which Vq can play an important role. 
We saw in Section ^ that the entries of the first row and column of R^M^ 1 R* D are smaller in 
the cases 1 and 2. If we allow cancellations between these entries of order A and the mixing 
parameter \J \[\, we might expect some qualitatively different results. In this more general case, 
we use fi', v 1 and p' to parameterize Vq rather than Lp. Including the effects of L e , we no longer 
have the symmetry constraints (/i, z/, p) ~ (//, z/, p'). 
We consider case 1: 



R^M^R 



D 



/ 



Here, although the unitary rotations remain relatively close to unity, the rotation parameter 
p! \f \ may be large enough to cancel the experimental term A. Such cancellation is only possible 
if p! ~ 1.Q Under geometrical constraints, the (11) entry will be \ y while the (12) entry could be 



11 Since [i! includes contributions from L e , its coefficient \J r]/e should be \frn € ~Jm(( l if the charged lepton ratio 
is larger. However, since ^m e /m^ ~ O.f ~ \J rj/e under our assumptions, we keep our familiar notation. 
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much smaller than — . One can proceed to analyze the mixing parameters in R D as in Section £2[ 
Due to the relatively large (11) entry, one finds that the constraints 

are required to preserve the small (12) and (13) entries. Since this requires \i -C y! ~ 1, the 
charged lepton rotations would have to be significantly larger than those from the Dirac neutrino 
matrix, at least for the 1-2 mixing. This would suggest an approximately geometric structure in 
the charged lepton matrix and a Dirac neutrino matrix with very small first generation mixing. 
Thus, the Dirac neutrino matrix would have a more restricted form than the cascade hierarchy 
derived for the simpler case without cancellations. Unless some additional information prompts 
us to favor these textures for M v , M D and M e , there is no compelling reason to further pursue 
this route. 

In the second case, where (M v )u ~ 0, we find that the constraints are the same as those 
listed in Eq. (|25|), i.e., the same as we found for this case without allowing for cancellations. 
These results hold because, regardless of how small A + p! ' \p\ may be, we retain the same 
relative hierarchy between the first generation entries and the same hierarchy in the 2-3 block, 
cf. Eq. dH). 

We conclude that these potential cancellations have little effect on our previous considerations. 

10 Outlook 

Barring cancellations or additional flavor symmetries, the observed pattern of neutrino mass 
splittings and mixing angles leads us to two related propositions for simple model building in 
the general context of a grand-unified theory with type-I seesaw mechanism. The first is a 
double hierarchy, with respect to the hierarchy of the Dirac matrix, Mo, in the effective heavy 
neutrino matrix Mr. The second, contingent upon the first, is a cascade structure in Mp, or a 
texture which is even closer to diagonal. These conclusions follow only from the structure of the 
type-I seesaw formula, together with the observation that the experimental neutrino data most 
naturally arise from an approximately democratic effective light neutrino matrix. If the neutrino 
masses obey a normal hierarchy, i.e., mi < m 2 ~ ^/Am^~ l <C m 3 ~ a/ Am^ tm , it is possible to 
relax these constraints, but it remains true that Mr should have an enhanced hierarchy and Mo 
should have a sub-geometrical structure. Moreover, in this case some approximate symmetry 
must exist to generate a second large mixing angle and a hierarchy consistent with experiment. 

These conclusions are rather general and not restricted to the specific model with small 
representations outlined in Section Q They hold for hierarchical, symmetric matrices, up to 
factors of order one. In light of the quark and charged lepton mass hierarchies, it is natural for 
Mo to be hierarchical. In particular, this matrix is closely related to the up quark matrix in 
many GUT models. Family symmetries will also tend to engender such relations. In Section [7| 
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we showed that even adding an operator which ostensibly only contributes to the Dirac neutrino 
matrix does not necessarily relax our conclusions. 

Can we implement these textures in a complete model? We discussed a scenario with a 
£7(1) x ^2 x Z2 flavor symmetry, where we generated a cascade structure for the Dirac matrices 
through the Froggatt-Nielsen mechanism. A double hierarchy in Mr is natural if it is an effective 
operator generated by integrating out singlets coupled to 16jl6#, where this coupling has an 
eigenvalue hierarchy similar to that in Mn (cf. Ref. [22]). However, this structure led to problems 
in the quark sector. We have seen that the relatively large Cabibbo angle implies that the down 
quark matrix is not purely cascade-like, although a cascade structure in the third generation is 
supported. This does not necessarily conflict with a fully cascade pattern in M D , but it requires 
a somewhat more complicated picture than the simple model described above. Furthermore, 
in our specific model, we rely on the antisymmetric operator 16jl6jl0#45# to differentiate the 
down quark and charged lepton matrices. This implies that its contributions cannot be too small. 
Since it also contributes to the up-quark and neutrino matrices, it becomes difficult to reconcile 
a cascade structure in these matrices with the strong up-quark hierarchy in a natural way. 

Lopsided models may provide us a way out of these potential difficulties. Compared with 
a cascade pattern, they are equally compatible with the CKM matrix. For the purposes of 
mass fitting, lopsidedness slightly relaxes the need for large off-diagonal contributions from 
16jl6jl0#45#. More importantly, a lopsided charged lepton matrix introduces large rotations 
which contribute to the PMNS matrix. If these are primarily responsible for one or both of the 
large mixing angles, it is possible to reduce the pull towards a double hierarchy in Mr. This in 
turn can relax the constraints that lead us to a cascade structure for Mb and so for M u . Exactly 
how much lopsidedness can obviate the need for a double hierarchy remains an open question. 
The atmospheric mass splitting remains small compared to the quark mass splittings, irrespec- 
tive of the origin of the large mixing angles. This will tend to require an enhanced hierarchy in 
at least part of Mr. Additionally, while it is technically possible that most or all of the PMNS 
structure comes from charged lepton rotations, we must ask how much can be done in a natural 
way. For example, as discussed at the end of Section [|, a small value for # 13 precludes generically 
large mixing from lopsidedness in all generations. 

This brings us to the nature and origin of #13 in general, which we have not addressed in detail 
in this paper. We chose to leave this an open question in light of the current uncertainty in the 
size of 8i 3 : only an upper bound is known. While it is clear that the solar and atmospheric mixing 
angles are large compared with those in the quark sector, #13 may or may not be comparatively 
small. Actually, the experimental upper limit, approximately 10°, is of the same order as the 
Cabibbo angle. This is large enough that its smallness compared to the other neutrino angles 
may be explained by normal fluctuations of order one parameters without violating our sense of 
naturalness ||. However, if #13 is significantly closer to zero we should seek some more robust 
explanation. For the forms of M v listed in Section ^|, this would require a symmetry closely 
relating various matrix elements. Another possibility arises for partially lopsided matrices: if 
one large mixing angle arises from the charged lepton sector and the other from M v then it is 
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natural to preserve a small third angle. Clearly, it is important to determine the order of #13. 

In summary, a combination of partially lopsided and partially cascade matrices, in conjunction 
with an enhanced hierarchy in Mr, seems to be the most natural route to explain the generic 
features of the quark and lepton data in a grand-unified model. The details of a complete model 
remain to be worked out, but our conclusions follow from a fairly general framework. It will be 
interesting to see if a workable model can be obtained with relatively simple family symmetries 
and what consequences there might be for experimental predictions. 
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A Derivation of Expanded Seesaw Formula 

In this appendix, we derive the extended seesaw formula, given in Eq. fl5E|). As mentioned in 
Section [7|, it is crucial that terms quadratic in M M do not appear. The formula was originally 
derived in Ref. through a slightly different calculation. 
As displayed in Eq. we propose the operators 

W a = Myl6jT6ifSj + M^.ieJe^S,- + m s (Ms)^ SSj ■ 

We integrate out the singlet fields, S, by taking a partial derivative and setting it equal to zero, 



dW s 
~dS~ 



0: Si 



2 m s 



MijWiWH + M' ij 16 t 16' H {Ms 1 ). 



(64) 



Plugging this into our initial equation yields 



W, 



cS 



Am. 



16, 



MMg X M 



T 



1Q H 1Q H + M M* X M 



T'T 



+ 2 [M Mg M 16^16^ 



16, 



(65) 



Now we let the Higgs fields acquire their GUT and weak scale vevs and we include the Dirac 
term uM^N, where v and N are the left and right-handed neutrinos, respectively. Suppressing 
the generation indices, we obtain 

1 



W, 



N 



4m s 



j 2 N ( MM^M T ) N + u' 2 u( m'Ms 1 M iT ) u + 2u'vu (m'm^M 7 ) N 



+ uvM D N , 

We extremize with respect to iV and find 
dW N 



dN 



= : N 



— 7T M D [M 



MqM 



MM 



v . 



(66) 



(67) 
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Inserting this into the last equation and performing a little algebra gives the amended seesaw 
formula 



Wf ~ |m d (M~ 1 ) T ' M s M~ l M T D - 1 U ' V 



2 um s 



M'M 1 Ml + M D (m'm lX ' 



u 2 m 



s 



V 2 



(68) 
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